AD-A144  630  ANALYTIC  MODELS  FOR  RADIATION  INDUCED  LOSS  IN  OPTICAL  1 
FIBERS  I<U)  ROME  AIR  DEVELOPMENT  CENTER  GRIFFISS  AFB  NV 
J  A  WALL  APR  84  RADC-TR-84-71 

UNCLASSIFIED  F/G  18/6  NL 


AD- A 144  630 


RADC-TR-84-71 
In-House  Report 
April  1984 


ANALYTIC  MODELS  FOR  RADIATION 
INDUCED  LOSS  IN  OPTICAL  FIBERS  I 


James  A.  Wall 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


>1  ^ 


\  AUC  2  t  '-34  .  ;{, 


A 


ROME  AIR  DEVELOPMENT  CENTER 
Air  Force  Systems  Command 
Griffiss  Air  Force  Basef  NY  13441 


84  08  20  16  4 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Office  (PA)  and 
is  releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS 
it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

RADC-TR-84-71  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 

i 


BOBBY  L.  BUCHANAN,  Chief 

Radiation  Hardened  Electronics  Technology  Branch 
Solid  State  Sciences  Division 


APPROVED: 


/ 


HAROLD  ROTH,  Director 
Solid  State  Sciences  Division 


FOR  THE  COMMANDER: 


JOHN  A.  RITZ 

Acting  Chief,  Plans  Office 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC  mailing 
list,  or  if  the  addressee  is  no  longer  employed  by  your  organization,  please 
notify  RADC  (  F.SR  )  Hanscom  APB  MA  01731.  This  will  assist  us  in  maintaining  a 
current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices 
on  a  specific  document  requires  that  it  be  returned. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


j  REPORT  DOCUMENTATION  PAGE  j 

1«.  REPORT  SECURITY  CLASSIFICA1  ION 

Unclassified 

lb  RESTRICTIVE  MARKINGS  | 

N/A  | 

2a  SECURITY  CLASSIFICATION  AUTHORITY 

N/A 

3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

2b.  OE CLASS! F  1  CATION /DOWNGRADING  SCHEDULE 

LNJA 

U.  PERFORMING  ORGANIZATION  REPORT  NUMBERS) 

5.  MONITORING  ORGANIZATION  REPORT  NOMBER(S) 

RADC-TR -84-71 

6a  NAME  OF  PERFORMING  ORGANIZATION 

7a.  NAME  OF  MONITORING  ORGANIZATION 

Rome  Air  Development  Center 

Rome  Air  Development  Center  (ESR)  1 

6c.  AOORESS  (City.  Stott  and  /.IP  Code/ 

7b  AOORESS  (City.  State  and  /IP  Code l 

Hanscom  AFB 

Massachusetts  01731 

Hanscom  AFB 
Massachusetts  01731 

8a  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 

9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER  | 

Rome  Air  Development  Center 

_ 

8c  ADDRESS  <Cit\.  State  and  /.IP  Code/ 

10  SOURCE  OF  FUNDING  NOS.  1 

Hanscom  AFB 

Massachusetts  01731 

PROGRAM 

element  no 

61102 

PROJECT 

NO 

2306 

TASK 

NO 

J3 

[11  TITLE  'Include  A*curifv  Classification) 

[ANALYTIC  MODELS  FOR  RADIATION  Icontd) 

12  PERSONAL  AUTHORISI 


.James  A.  Wall 


IU»i!HWrtOTMMi 

■j 

15  PAGE  COUNT  ft 

_ 

16  SUPPLEMENTARY  NOTATION 


17 

COSATl 

COOES 

Jfl  SUBJECT  TERMS  >t 'on  tm  up  on  reverie  if  necestar^  and  identify  by  block  number 1 

F  lELD 

GROUP 

SUB  GR 

Optical  fibers 

20 

12 

Radiation  effects 

Analytic  models 

19  ABSTRACT  Con  l  in  up  on  rrieru-  if  necenars  and  identify  6\  blot  *  number: 


Analytic  models  to  describe  the  radiation  induced  loss  in  optical  fibers  as  a  function  of 
dose  were  sought.  A  simple  power-law  function  was  tried,  and  it  w'as  found  to  give  excel¬ 
lent  to  good  fits  to  data  on  2  5  fibers  having  a  wide  range  of  radiation  sensitivities. 

Using  this  function,  only  2  empirical  constants  are  needed  to  calculate  the  radiation 
induced  loss  in  an  optical  fiber.  This  is  useful  in  categorizing  the  radiation  hardness 
of  optical  fibers  and  in  systems  d esien. 

N 


20  DISTRIBUTION  AVAILABILITY  QF  ABSTRACT 

21  ABSTRACT  SECURITY  CLASSIFICATION  1 

UNClASSiF'EO  UN  l  1  Mi  TEQ  —  SAME  AS  APT  _  DTtC  USERS  H." 

Unclassified 

22a  NAME  OF  RESPONSIBLE  iNOiViDuAl 

221.  TELEPHONE  NUMBER 

Include  1  Yfl  Code 

22c  OF  F  iCE  SYMBOL 

James  A.  Wall 

(617)  861-4031 

RADC/ESR 

|  James  A.  W  all 

DO  FORM  1473.  83  APR 


EDiTiQN  OF  1  jAN  73  iS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


Contents 


1.  INTRODUCTION  5 

2.  PROCEDURE  6 

3.  RESULTS  8 

4.  DISCUSSION  10 

REFERENCES  12 


Illustrations 


1.  Example  of  Excellent  Fits  of  Eq.  (1)  to  Fiber  Optic  Radiation 

Induced  Loss  Data  9 

2.  Example  of  "Good"  Fits  of  Eq.  (1)  to  Fiber  Optic  Radiation  Induced 

Loss  Data  9 


■  List 


Analytic  Models  for  Radiation  Induced 
Loss  in  Optical  Fibers  I 


I.  INTRODUCTION 


During  a  program  on  the  development  of  radiation-hardened  optical  fi- 
12  3 

bers,  *  *  a  large  amount  of  data  was  accumulated  on  the  radiation  induced  loss 

in  optical  fibers  as  a  function  of  dose.  Data  taken  during  the  initial  phases  of  the 
program  showed  a  linear  increase  in  induced  loss  with  increasing  dose  so  that  the 
radiation  hardness  of  the  fibers  could  be  characterized  by  the  slopes  of  plots  of  the 
data,  that  is,  dB/km-kilorad.  However,  as  the  accuracy  of  the  radiation  test  pro¬ 
cedure  was  improved,  significant  departures  from  linearity  were  observed  in  the 
data,  particularly  at  doses  below  about  5  kilorads.  Nonlinear  behavior  of  induced 

loss  in  optical  fibers  as  a  function  of  dose  had  been  reported  by  others  (see  Frie- 
4 

bele,  for  example).  However,  their  fiber  irradiations  were  performed  at  much 
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higher  dose  rates  than  ours,  and  the  nonlinearities  they  observed  could  be  attri¬ 
buted  to  failure  to  establish  dynamic  equilibrium  between  the  trapping  (light  ab¬ 
sorbing)  centers  and  the  electron-hole  pairs  generated  during  irradiation.  The 
failure  to  establish  dynamic  equilibrium  during  the  irradiations  resulted  in  simul¬ 
taneous  generation  and  annealing  of  light  absorbing  centers  during  irradiation. 
These  results  could  account  for  nonlinearity  of  the  data.  This  could  be  observed 
as  partial  annealing  of  the  induced  loss  over  a  period  of  several  minutes  following 
the  irradiations,  a  feature  which  our  data  did  not  show. 

An  analytic  expression  to  describe  the  radiation  induced  loss  in  optical  fi¬ 
bers  as  a  function  of  dose  is  important  in  simplifying  the  categorization  of  the  ra¬ 
diation  hardness  of  fibers  as  well  as  providing  the  design  engineer  with  a  means 
of  determining  fiber  performance  without  reference  to  data  plots.  Since  the  linear 
description  of  radiation  response  could  no  longer  be  considered  a  reasonable  ap¬ 
proximation,  other  functional  relationships  were  sought.  It  was  found  that  some 
workers  had  reported  a  dependence  of  the  induced  loss  on  the  square  root  of  the 
dose'*’  Although  the  data  reported  by  these  workers  appeared  only  to  ap¬ 

proximate  the  square  root  of  dose  dependence,  it  was  decided  to  attempt  a  power- 
law  fit  to  our  data. 

2.  IMHM  KIH  UK 

The  data  on  radiation  induced  loss  vs  dose  of  25  optical  fibers  were  fit  to  the 
expression 

L.  =  ADn  (1) 

where  L  =  induced  loss  in  dB/km 
D  =  dose  in  kilorads 
A,n  -  constants 

A  least  squares  fitting  procedure  was  performed  using  a  hand  calculator. 
The  procedure  is  described  in  most  elementary  textbooks  on  statistics  and  is 
outlined  here  for  information. 

5.  Share,  S. ,  and  IVasilik,  J.  (1979)  Radiation  effects  in  doped-silica  optical 

waveguides,  IEEE  Trans,  on  Nucl,  Sci,  26:4802. 

6.  Gur'yanov,  A.N.,  et  al  (1978)  Radiation-optical  stability  of  low -loss  glass-fi¬ 

ber  waveguides,  Sov.  J.  Quantum  Electronics  8:1401. 

7.  Gur'yanov  et  al  (1979)  Radiation-optical  stability  of  low-loss  glass-fiber 

waveguides,  Spy,  J.  Quantum  Electronics  9,:768. 


Eq,  (1)  was  linearized  by  taking  the  logarithm  of  both  sides 

Log  L  =  Log  A  -  n  Log  D 

define  y  =  Log  L 
x  =  Log  D 
a  =  Log  A 

The  required  constants  are  then  given  by: 

xy  -~x  •  y 


2  2 
x  -  x 


a  =  y  -  n  •  x  ;  A  =  antilog  a  (4) 

where  "x  =  average  value  of  the  x's  (Log  D's) 
v"  =  average  value  of  the  y's  (Log  L's) 

xy  *  average  value  of  the  products  of  corresponding  x's  and  y's 
x^  =  average  value  of  the  squares  of  the  x's 

Eq.  (3)  is  most  frequently  presented  in  terms  of  summations  over  the  indivi¬ 
dual  data  points  in  textbooks,  but  the  form  given  here  is  more  practical  for  use 
with  basic  hand  calculators. 

Although  a  correlation  coefficient  can  be  calculated  as  a  measure  of  the 
goodness  of  fit"  for  each  data  set,  a  more  definitive  evaluation  of  how  well  the 
fitted  analytic  expression  agrees  with  the  data  is  obtained  by  calculating  the  root- 
mean-square  (rms)  of  the  deviations  of  values  calculated  using  the  fitted  expres¬ 
sion  from  corresponding  values  of  the  data.  We  define  the  rms  of  the  deviation  as 


uV'2 


.  .2  .  <Lmi-Lct»2 

where  A  =  L 


Lmi  =  induced  loss  measured  at  dose  i 
Lei  =  induced  loss  calculated  at  dose  i 
k  =  number  of  data  points 
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The  fits  of  the  data  on  25  optical  fibers  to  Eq.  (1)  ranged  from  excellent  to 
good.  To  illustrate  the  meanings  of  these  qualitative  fit  evaluations.  Figure  1 
shows  examples  of  two  excellent  fits,  one  example  for  a  relatively  radiation  sen¬ 
sitive  fiber  and  one  for  a  more  radiation  "hard"  fiber.  As  can  be  seen  from  the 
figure,  the  curves  calculated  from  the  fits  are  essentially  identical  with  the  plots 
of  the  actual  data  points.  Examples  of  the  poorer  fits,  which  we  can  still  call 
good,  are  shown  in  Figure  2  for  a  similar  pair  of  fibers.  Although  the  fit  for  the 
more  radiation  sensitive  fiber  tends  to  overestimate  the  induced  loss  at  higher 
doses,  and  the  fit  for  the  more  radiation  hard  fiber  underestimates  the  induced 
loss  at  the  higher  doses,  both  fits  are  quite  acceptable  for  the  lower  doses. 

For  a  more  quantitative  evaluation  of  the  quality  of  the  fits,  the  rms  deviation 
of  induced  losses  calculated  using  the  fitted  Eq.  (1)  from  the  measured  induced 
losses  were  calculated  using  Eq.  (5).  The  average  rms  deviation  for  all  the  fibers 
was  0.  698  dB/km  with  a  standard  deviation  of  0.484  dB/km.  Considering  that  the 
estimated  accuracy  of  the  data  was  about  +2  dB/km,  this  means  Eq.  (1)  can  repre¬ 
sent  the  induced  loss  as  a  function  of  dose  very  well.  In  fact,  for  17  of  the  25  fi¬ 
bers  considered,  the  rms  deviation  of  the  fits  from  the  data  was  less  than  1 
dB/km. 

The  quality  of  the  fits  is  surprising  considering  that  a  fairly  wide  range  of 

radiation  sensitivities  was  covered  by  the  optical  fibers  observed.  This  range  of 

12  3 

sensitivities  is  illustrated  by  Wall.  '  '  To  summarize,  the  group  consisted  or 

fibers  having  germanium-silicate  cores  doped  to  varying  degrees  with  boron, 
phosphorous,  antimony,  and/or  cerium  in  different  amounts  and  combinations 
giving  a  range  of  radiation  sensitivities  of  about  an  order  of  magnitude.  Of  the  two 
constants  required  to  fit  the  data  to  Eq.  (1),  A  was  the  most  variable.  For  all  the 
fibers,  the  average  value  of  A  was  6.99  with  a  standard  deviation  of  4.  54. 

The  higher  values  of  A  could  be  identified  with  fibers  having  phosphorous 
added  to  the  cores.  The  average  value  of  A  for  these  fibers  was  10.  32  with  a 
standard  deviation  of  4.  11.  For  the  fibers  with  boron  added  to  the  core,  the  aver¬ 
age  of  A  (omitting  one  unusually  radiation  sensitive  fiber  from  the  group)  was  4.23 
with  a  standard  deviation  of  I.  16.  In  both  cases,  the  variations  in  A  could  only  be 
roughly  correlated  with  the  addition  of  antimony  or  cerium  to  the  fiber  cores.  The 
values  obtained  for  n  showed  much  less  spread  than  those  for  A.  For  all  25  fibers, 
the  average  value  of  n  was  0.  860  with  a  standard  deviation  of  0.  073.  For  the  phos¬ 
phorous  doped  fibers,  the  average  value  of  n  was  0.923  with  a  standard  deviation 
of  0.045.  l  or  the  boron  doped  fibers,  the  average  value  was  0.  826  with  a  stand¬ 
ard  deviation  of  0.  074.  There  appeared  to  be  a  somewhat  better  correlation  be¬ 
tween  the  values  of  ti  and  the  addition  of  antimony  or  cerium  to  the  fiber  cores. 


Both  these  dopants  tended  to  lower  the  radiation  sensitivities  of  some  fibers,  and 
this  generally  corresponded  with  lower  values  of  n. 

The  foregoing  applies  to  fibers  irradiated  at  room  temperature,  A  limited 
amount  of  data  was  also  obtained  on  some  of  these  fibers  irradiated  at  tempera¬ 
tures  of  -55°  C  and +125°  C.  Eq.  (1)  could  also  be  fit  fairly  well  to  these  data. 
Because  of  the  small  amount  of  data,  not  much  can  be  said  about  these  results  ex¬ 
cept  that  the  same  analytic  expression  (with  different  values  for  the  constants)  ap¬ 
parently  can  be  applied  over  these  temperature  extremes. 


».  DISCISSION 

It  is  surprising  that  such  a  simple  expression  as  Eq.  (1)  can  be  used  to  des¬ 
cribe  the  radiation  induced  loss  as  a  function  of  dose  in  a  variety  of  optical  fibers 
of  different  radiation  sensitivities  and  compositions.  At  present,  we  cannot  think 
of  any  physical  reason  why  this  is  so. 

Those  who  have,  justifiably  or  not,  invoked  the  square  root  of  dose  depend¬ 
ence  of  induced  loss  that  prompted  this  study  have  given  two  reasons  why  a  simple 
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power  law  should  be  applicable  to  optical  fibers.  The  first  rationale  '  is  based 
on  a  relatively  simple  model  used  to  explain  the  creation  of  F  centers  (color  cen- 

g 

ters)  in  crystalline  lithium  fluoride  by  X  rays.  It  seems  unlikely  that  this  work 
could  be  simply  extended  to  the  amorphous  compound  glass  from  which  optical  fi¬ 
bers  are  made.  Of  course,  the  model  does  not  apply  to  optical  fibers  since  it  pre¬ 
dicts  a  strict  square  root  of  dose  dependence,  which  is  apparently  not  the  usual 
case  found  for  the  fibers. 

5 

The  second  rationale  is  based  on  the  assumption  that  annealing  of  the  induced 

Q 

loss  occurs  during  irradiation.  One  calculation  has  been  made'"  based  on  this  as¬ 
sumption  for  a  plastic  optical  fiber.  Transient  data  on  the  recovery  of  induced 
loss  following  an  X-ray  pulse  was  found  to  be  proportional  to  t_rn  where  t  is  the 
time  and  m  a  constant.  This  was  combined  with  an  assumed  induced  loss  genera¬ 
tion  function  to  obtain  a  power  law  dependence  on  time  and  dose-rate  that  agreed 
reasonably  well  with  the  induced  loss  data.  However,  this  agreement  was  probab¬ 
ly  fortuitous.  From  our  own  measurements  on  transient  radiation  effects  in  opti- 


8.  Durand,  P. ,  Farge,  V,,  and  Lambert,  M.  (1969)  The  creation  of  F  centers  in 

lithium  flouride  between  77°  and  600°  K  and  their  interpretation  by  a  re¬ 
combination  model  of  interstitial -vacancies,  J.  Phvs.  Chem.  Solids 
30:1353. 

9.  Kalma,  A.H.,  and  Hardwick,  W.H.  (1976)  Radiation  testing  of  a  fiber  optics 

data  link,  IEEE  Trans,  on  Xucl.  Sci.  23:1769. 
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cal  fibers  and  the  measurements  of  others,  it  is  evident  that,  in  general,  the 
recovery  of  induced  loss  is  not  a  simple  function  of  time.  This  is  probably  because 
of  the  presence  of  different  trapping  centers  with  different  trap  release  times  in  a 
given  fiber.  Also,  since  our  fibers,  irradiated  at  much  lower  dose-rates,  showed 
no  recovery  following  irradiation,  the  "recovery  during  irradiation"  explanation 
is  not  considered  applicable. 

In  spite  of  the  lack  of  physical  justification  for  the  power  law  dependence  of 
radiation  induced  loss  on  dose  in  optical  fibers,  the  possibility  of  describing  the 
induced  loss  data  in  terms  of  only  two  constants  is  a  worthwhile  finding  for  compi¬ 
ling  data  on  the  radiation  hardness  of  optical  fibers  and  in  opto-electronics  system 
design.  However,  as  with  any  such  simple  analytic  model,  caution  must  be  used 
in  its  application  because  the  values  for  the  required  constants  obtained  can  be  ex¬ 
pected  to  be  dependent  on  irradiation  conditions  such  as  dose-rate  and,  as  we 
already  know,  temperature. 


10.  Wall,  J.A.,  Posen,  H.,  and  Jaeger,  R.  (1981)  Temperature  Response  of 

germanium  phosphosilicate  optical  fibers  under  irradiation,  in  Physics  of 
Fiber  Optics,  Advances  in  Ceramics,  B.  Bendow  and  S.  S.  Mitra,  Eds., 
American  Ceramic  Society  2:398. 
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5  MISSION  5 

i  of  * 

^  Rome  Air  Development  Center  ? 

n  MOC  plan*  and  execute*  research,  de.vzi.opme.nt,  test  and  V 

•  selected  acquisition  pn.ogn.am  in  support  o (J  Command,  Con#Lo£  O’ 
5  C ormunicatio n*  and  Intelligence  (C3 1)  activities* .  Technical  ? 

•  and  engineering  4uppo4t  within  area. s  0f{  tecfintcai.  competence  C 

S  -C4  provided  to  ESV  Pn.ogn.am  Voices  (POs)  and  other  ESV  § 

^  element a.  The  pnincipal  technical  mission  area*  one  \ 

S  communication* ,  electromagnetic  guidance  and  control,  sun- 

D  vdliance  ground  and  aeAo^pace  object*,  intelligence  data  \ 

•  cottectcon  and  ftandCtng,  tn^oomatxon  4  intern  tec/inotogg,  0 

k  ionospheric  propagation,  solid  a  tote  -science* ,  microwave 
t  physic*  and  electronic  reliability,  maintainability  and 
S  compatibility. 


